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The properties of matter at ultra-high densities, low temperatures, and with a significant asym- 
metry between protons and neutrons can be studied exclusively through astrophysical observations 
of neutron stars. We show that measurements of the masses and radii of neutron stars can lead to 
tight constraints on the pressure of matter at three fiducial densities, from 1.85 to 7.4 times the den- 
sity of nuclear saturation, in a manner that is largely model-independent and that captures the key 
characteristics of the equation of state. We demonstrate that observations with 10% uncertainties of 
at least three neutron stars can lead to measurements of the pressure at these fiducial densities with 
an accuracy of 0.11 dex or ~ 30%. Observations of three neutron stars with 5% uncertainties are 
sufficient to distinguish at a better than 3a confidence level between currently proposed equations 
of state. In the electromagnetic spectrum, such accurate measurements will become possible for 
weakly-magnetic neutron stars during thermonuclear flashes and in quiescence with future missions 
such as the International X-ray Observatory (IXO). 

PACS numbers: 97.60. Jd, 26.60.Kp, 21.65.-f 



I. INTRODUCTION 

Matter in the deep interiors of neutron stars is char- 
acterized by low temperatures, large chemical potentials, 
and a significant asymmetry between the number of neu- 
tron and protons. In these conditions, several phenomena 
may take place, such as the creation of stable bosons and 
hyperons or a phase transition to quark matter. In ad- 
dition, at these high densities, the symmetry energy of 
nucleonic matter plays a dominant role in determining 
its microscopic properties. The large uncertainty in un- 
derstanding these phenomena has led to a wide range of 
models for the equation of state of neutron-star matter, 
which cannot be distinguished by hadron collision exper- 
iments or with observations of the early universe [l|, 0] . 

The relation between the pressure P and density p of 
neutron-star matter determines the macroscopic proper- 
ties of the stars and, in particular, their radii R, masses 
M, and moments of inertia I. In fact, there is a unique 
map between the microscopic P — p relation and the 
macroscopic M — R one [3]. Astrophysical observations 
of neutron stars aim to exploit this mapping and invert 
it in order to measure their equation of state. 

The mass and radius of a neutron star of a given cen- 
tral density, pc, depends on the entire relation between 
pressure and density up to pc, because of the density 
gradient from the center of the star to its surface. If 
the mass-radius relation could be traced observationally 
with a large number of neutron stars spanning the entire 
range of masses (e.g., 0.2-2 Mq), then the mapping could 
be formally inverted to generate the complete equation 
of state 3. However, the astrophysical formation chan- 
nels for neutron stars limit them to masses larger than 
~ I.2M0, rendering this approach impractical. 

Prakash and Lattimer [3| studied extensively the prop- 
erties of neutron-star models for a large range of equa- 
tions of state and concluded that the predicted radii 
depend primarily on the pressure at a fiducial density, 
comparable to the nuclear saturation density of pns ~ 



2.7 X 10^"* g cm More recently. Read et al. [5] showed 
that not only the approximate radii of neutron stars, but 
a large number of observables depend on a small number 
of parameters that characterize the equation of state. 

In this paper, we show that the masses and radii of neu- 
tron stars act as tracers of the pressure of neutron-star 
matter at three fiducial densities. In other words, the 
complete mass-radius relation to high accuracy can be 
reproduced for all proposed equations of state, when the 
pressure at these three densities is specified. Therefore, 
a one-to-one mapping exists between the pressure at the 
three densities and the mass-radius relation, which can be 
inverted even if the observed neutron stars are confined 
to a small range of masses. We present the formalism for 
this inversion that takes into account the uncertainties 
in the mass and radius measurements and show that ob- 
servations of three neutron stars with a 10% uncertainty 
can lead to a measurement of the pressure of ultra-dense 
matter at these densities with comparable accuracy. 



II. A MINIMAL REPRESENTATION OF 
NEUTRON-STAR EQUATIONS OF STATE 

Inferring the pressure of neutron-star matter at a large 
number of densities from astrophysical measurements 
would ideally generate the equation of state in a model 
independent way. However, the range of core densities of 
neutron stars with masses > I.2M0 is in fact not very 
large and, therefore, the pressure at a small number of 
fiducial densities pi governs their observable properties. 
Choosing two fiducial densities, setting them to the fixed 
values pi = l.SSpns and p2 = 2pi, and assuming a piece- 
wise polytropic relation between them has been shown to 
accurately reproduce the complete pressure-density rela- 
tion for a large sample of proposed equations of state [1] . 

We explored the predicted mass-radius relations for the 
parametric equation of state of Ref. [3] in which the poly- 
tropic indices between the fiducial densities are left as 
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free parameters. We found that different combinations 
of polytropic indices lead to practically indistinguishable 
mass-radius relations for the neutron star models. Con- 
sequently, this choice leads to a significant degeneracy 
between the inferred uncertainties of the various model 
parameters when aiming to reconstruct the equation of 
state from a set of mass-radius measurements (see §111) . 
This degeneracy is greatly reduced if we choose, instead, 
the values of P{pi) and P{p2) to be the free parameters 
and describe the equation of state at densities larger than 
P2 by a polytrope with normalization set by the pressure 
P{p3) at a third fiducial density pz > P2- Because the 
density /03 only serves to specify the polytropic index at 
densities larger than p2, its exact value does not affect 
the final result. Hereafter, we set r/ios = 2 x p2. Fol- 
lowing Ref we supplement this with the equation of 
state SLy Q for the outer layers of the neutron star for 
densities below a pQ, which we leave as a free parameter. 

Specifically, we use hereafter the following parametric 
equation of state of neutron-star matter, which depends 
on four parameters po, Pi = P{pi), P2 = P{p2), and 
P:i=P{p:i). 

(i) For p < po, we interpolate between the tabulated val- 
ues of the SLy equation of state [Sj] . We define Pq and cq 
to be the pressure and energy density at po, respectively. 
(a) For po < p < pi, we define 



Fi 



log(Pi/Po) 



log(pi/po) 

so that the pressure in this density range is given by 



P^Pi[ — 

.Pi 



and the energy density is 

e = {l + ai)p + 

where 
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(Hi) For pi < p < P2, we set 

_ log(P2/Pl) 
log(p2/Pl) 

SO that the pressure in this density range is given by 
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and the energy density is 

e== (l + a2)p + 
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FIG. 1: Mass-radius relations for neutron stars, computed us- 
ing the complete (solid lines) as well as the parametric (dashed 
lines) forms of three equations of state. 



TABLE 1: Best-fit EOS parameters for three sample equations 
of state. 



EOS 


log Po 


log Pi 


log Pa 


logPa 


EPS 


14.30 


34.283 


35.142 


35.925 


GSl 


14.85 


34.504 


34.884 


35.613 


APS 


14.30 


34.392 


35.464 


36.452 



where 



02 = ai 



Pi 



Pi 



(Fi-i)pi (F2-i)pr 

(iv) For p > p2, we set 

_ l0g(P3/P2) 
^ I0g(p3/p2) 

SO that the pressure in this density range is given by 



(8) 
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P = P2 I ^ 

P2 



and the energy density is 
e = (H-a3)p + 



Fs - 1 \P2 



where 

03 = 02 + 



Pi 



P2 



P2 



(F2 - l)p2 \Pl J (F3 - l)p2 ' 



(10) 



(11) 



(12) 



Figure [T] shows the mass-radius relations obtained by 
integrating the Tolman-Oppenheimer-Volkoff equations 
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FIG. 2: The dependence of the calculated mass-radius relations of neutron stars on the parametric equation of state. The first 
three panels show the change in the predicted relation when the values of the parameters Pi, P2, and P3 are varied by 25% 
in each direction away from the best-fit values for the equation of state FPS. The last panel shows the change in he predicted 
relation when the value of the parameter po is reduced to 1/2 and 1/4 of its best-fit value for the same equation of state. 



with the complete (solid lines) as well as the parametric 
(dashed lines) forms of a sample of proposed equations 
of state, chosen to represent a wide range of physical 
conditions (see Table 1 for the parameters). The three 
examples also span a large range in neutron star masses 
and radii. For neutron stars with masses > IMq, the 
deviation in radius between the models with the complete 
and parametric equations of state is always < 2%. The 
same is true for all the equations of state considered in 
Ref. i. 

In order to get an intuitive understanding of the ef- 
fect of these parameters, we show in Figure [2] how the 
calculated mass-radius relations change when we vary 
each parameter in succession away from the best-fit val- 



ues for equation of state FPS. The value of the pressure 
at pi determines predominantly the characteristic radii 
of neutron stars; this is the effect described in Ref. [Ij. 
The value of the pressure at p2 regulates the slope of the 
mass-radius relation, at larger masses. The pressure at 
Pa sets the maximum mass of neutron stars. Finally, the 
value of the density poj where the parametric branch of 
the equation of state begins, only affects the mass-radius 
relation at small neutron-star masses. In the first three 
cases, the parameters are varied by 25% away from their 
best fit value in each direction while, in the last case, the 
density po is reduced to 1/2 and 1/4 of its original value 
in order to give rise to a visible difference in the resulting 
mass-radius relation. Because of the much weaker depen- 
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FIG. 3: Simulated values of the masses and radii of three 
neutron stars that obey the FPS equation of state (solid dots). 
The hatch-filled regions and the ellipses indicate 5% and 10% 
uncertainties in each measurement, respectively. 



dence of the observed properties of neutron stars on po, 
especially at the mass range > 1.2Mq of interest, we fix 
hereafter its value to log po — 14.30. 

The distinct dependences of the mass-radius relations 
on the values of the three parameters Pi, P2, and P3 
shown in Figure [2] will allow a determination of each pa- 
rameter from observations with minimal correlated un- 
certainties. Moreover, the rather strong dependence of 
the mass-radius relations on the parameters of the equa- 
tion of state guarantees that these parameters can be ac- 
curately measured from astrophysical observations. We 
discuss the expected uncertainties in the following sec- 
tion. 



III. EQUATION OF STATE PARAMETER 
ESTIMATION FROM MEASUREMENTS OF 
MASSES AND RADII 

The mass-radius relation of neutron stars that can be 
probed by astrophysical observations is determined to a 
large extent by the three parameters Pi, P2, and P3 of 
the parametric equation of state that we discussed in the 
previous section. This is because the value of the fourth 
parameter, po, only affects significantly the properties of 
neutron stars with masses that are too low to be astro- 
physically relevant. 

For a given equation of state, the particular realiza- 
tion of a neutron star is determined by its central density 
Pc- As a result, independent mass and radius measure- 
ments of N neutron stars require N + 3 parameters to 
be completely modeled, under the assumption that all 
stars follow the same equation of state. In other words. 
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FIG. 4: The probability density over the pressure P2 at 3.7 
times the nuclear saturation density, evaluated when the other 
two parameters of the equation of state are set to the best- 
fit values for equation of state FPS. This is the probability 
that the parametric equation of state reproduces the set of 
simulated data shown in Figure [3] that have 10% measure- 
ment uncertainties. The vertical dashed line indicates the 
true pressure at density p2 of the FPS equation of state. 



observations of N neutron stars yield 27V measurements, 
while modeling their properties requires A'^ -I- 3 parame- 
ters, leaving 2N — {N + 3) — N ~ 3 degrees of freedom. 
Thus, in order to obtain a unique solution, at least three 
neutron-star mass and radius measurements are required. 
Hereafter, we discuss the situation with iV = 3. 

Our aim is to invert the one-to-one map be- 
tween the parameter space of the observables 
{Ri,Mi,R2,M2,R3,Ms) to the one of the physi- 
cal quantities (Pi, P2, P3, Pc,i: Pc,2, Pc,3), where the 
indices on the masses, radii, and central densities 
refer to the three observed neutron stars. Because 
of the uncertainties inherent to the measurements, 
we actually need to convert a probability distribu- 
tion density over the parameter space of observables 
'Pohs{Ri, Ml, R2, M2, R3, M3) to one over the physical 
parameters P(Pi, P2, P3, Pc,i, Pc,2, Pc,3)- The particular 
values of the central densities for the three stars are not 
of particular interest and obtaining the final probability 
density over the equation-of-state parameters Pi, P2, 
and P3 would require marginalizing over the three central 
densities. In order to avoid this additional complication, 
we choose instead to transform the probability density 
from the parameter space of observables to one over 
(Pi, P2, P3, Ml, M2, M3) and then marginalize over the 
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FIG. 5: Contours of 68.2% (la) and 95.4% {2a) confidence levels for the values of pairs of parameters of the equation of state 
necessary to describe the simulated data shown in Figure [S] In each panel, the probability density is marginalized over the 
parameter not shown. The parameters Pi, P2, and Pg are the pressures at 1.85, 3.70, and 7.40 times the nuclear saturation 
density, respectively. In all panels, the corresponding true pressures for a sample of proposed equations of state is shown for 
comparison. The left panels correspond to a 5% uncertainty in the measurements, whereas the right panels correspond to a 
10% uncertainty. A 5% measurement error is sufficient to distinguish to a better than 3a level not only between equations of 
state that have different underlying physics but also between ones with similar mass-radius predictions. 



three neutron star masses. Formally, we write 

ViPi,P2,P3,A'h,M2,M3) = 

7'obs(i?l,Mi,i?2,M2,i?3,M3) 



Pl,P2,P3 



, (13) 



where the last term is the Jacobian of the transformation 
between the three observed radii and the three parame- 
ters of the equation of state. 

For each set of parameters of the equation of state 
Pi, P2 J -Pa I we obtain the relation between neutron-star 
radius and mass, R = R{M; Pi, P2, P3), by integrating 
the Tolman-Oppenheimer-Volkoff equation and use it to 



calculate numerically the Jacobian as 

Ri,R2,R3\ , /„ N , / dR 



J 



Pl,P2,P3 



det (Jij ) = det 



dP, 



(14) 



M, 



We then integrate equation (IT51) over the three neutron 
star masses to obtain the final probability distribution 
over the parameters of the equation of state 



P(Pi,P2,P3) = J j j dMidM2dM3 

P(Pl,P2,P3,Mi,M2,M3) . (15) 

In order to study the accuracy with which the equa- 
tion of state parameters can be inferred from astrophysi- 
cal measurements and to assess the potential correlations 
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between them, we simulated mass-radius data for three 
neutron stars that obey the FPS equation of state, cho- 
sen here as an example. In particular, we set Mi = 
1.50Mq, i?i ^ 10.7 km, M2 = 1.61Mq, R2 = 10.5 km, 
M3 = I.7OA/0, i?i = 10.2 km, and assigned to each 
measurement an independent Gaussian uncertainty with 
standard deviation equal to 5% and 10% of the corre- 
sponding central values. We show the simulated data in 
Figure |3] together with their uncertainties. Note that this 
is a very conservative and narrow range of measurements 
for neutron-star masses, as observed masses have been 
reported ranging from ~ 1.2Mq in radio pulsars Q to 
~ 2.IM0 in accreting neutron stars 0- We will discuss 
at the end of this section the improved constraints on 
the equation of state that can be obtained when a wider 
range of masses is considered. 

Using this set of simulated measurements, we inverted 
the mapping between neutron star mass-radius and the 
equation-of-state pressures at three fiducial densities and 
calculated the probability distribution over the three pa- 
rameters Pi, P2, and P3 according to equations (fT3|) to 
([T5|) . Because the resulting probability density is defined 
in a three dimensional parameter space, we show, in the 
following set of figures, projections and marginalizations 
of this function along different directions. 

In Figure |4l we show the normalized probability den- 
sity over the pressure P2 at the fiducial density p2, eval- 
uated when the other two parameters of the equation of 
state are set to the best-fit values for equation of state 
FPS (see Table |T| . The true value of the pressure P2 for 
equation of state FPS, shown by the vertical dashed line, 
is within 0.02 dex or 4.7% of the peak of the probabil- 
ity distribution. Moreover, the parameter P2 is tightly 
constrained with an la uncertainty of 0.11 dex or 29%. 

In Figure [5l we show the normalized two-dimensional 
probability distributions over different pairs of the equa- 
tion of state parameters, marginalized over the remaining 
parameter. In each panel, the two contours correspond 
to the 68.2% (Ict) and 95.4% (2cr) confidence levels. In 
all panels, the true pressures for a set of representative 
equations of state are shown for comparison. Finally, the 
left and right panels show the probability densities ob- 
tained for the assumed 5% and 10% uncertainties in the 
simulated data, respectively. 

The results shown in the same figure provide addi- 
tional justification for the use of three (as opposed to 
fewer) fiducial densities at which to measure the pres- 
sure of matter. In Ref. which showed that many pro- 
posed equations of state can be described by a piecewise 
polytropic model, this was necessitated by that fact that 
the potential production of bosons or hyperons at high 
densities or a phase transition to quark matter causes 
abrupt changes to the equation of state at multiple den- 
sities. Here, the use of three fiducial densities allows us to 
distinguish between equations of state, without requiring 
unattainable accuracy in the measurements. In the ex- 
amples shown in Figure [Sj while both equations of state 
FPS and GSl are consistent with the simulated data in 
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FIG. 6: The tightening of the constraints on the equation of 
state parameter Pi compared to the upper-left panel in Fig- 
ureO when a 1.3 Mq neutron star is included in the simulated 
data. 



the Pi — P2 parameter space, equation of state GSl can 
be excluded at a confidence better than 2 — Scr in the 
P2 — P3 space. 

Of the three pressures Pi, P2, and P3, the one at 3.7 
times nuclear saturation density, i.e., P2, is the most 
tightly constrained by the data. This is because of the 
particular range of masses that comprise our simulated 
data, for which the central density is comparable to p2- 
The pressure Pi at 1.85 times nuclear saturation den- 
sity determines mostly the equation of state outside of 
the cores of these neutrons stars and is, therefore, more 
tightly constrained from above (see the top panels of Fig- 
ure [5]). On the other hand, the pressure P3 at 7.4 times 
nuclear saturation density primarily affects the maximum 
allowed mass for neutron stars and is, therefore, more 
tightly constrained from below (see the bottom panels of 
Figure [5]). 

The above discussion suggests that allowing for a wider 
range of masses in the simulated data would have im- 
proved the constraints on Pi, if the range extended to 
lower neutron-star masses, and on P3, if the range in- 
cluded higher masses. We illustrate this for the case of 
Pi in Figure [6l where we have assumed that the three 
neutron stars have Mi = 1.31Mq, Ri — 10.95 km, 
M2 = I.5OM0, i?2 = 10.7 km, and M3 = 1.70Mq, 
Ri = 10.2 km and allowed for 5% measurement uncer- 
tainties. Comparison with the top- left panel of Figure [S] 
shows that the allowed range of values for the parameter 
Pi has shrunk significantly when a lower-mass neutron 
star is included. 

It is evident from Figures [5] and [6] that the uncertain- 
ties in the three pressures Pi, P2, and P3 are weakly 
correlated. In fact, there is no discernible correlation 
between P3 and the other two parameters. This moti- 
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FIG. 7: The constraints on the polytropic indices Fi and F2, 
which appear in the alternative parametrization of the equa- 
tion of state, obtained for the same simulated data as in Fig- 
ure |6] This parametrization leads to significantly correlated 
uncertainties even for the wide range of neutron star masses 
considered here and for 5% errors in the measurements. 



vated our choice of the three pressures as the parameters 
of the minimal representation of the equation of state 
as opposed to the polytropic indices between the corre- 
sponding densities (cf. Q). Using the parametrization 
with polytropic indices would have resulted in correlated 
uncertainties, such as those shown in Figure [7] for the 
case of Fi and F2, even for the most constraining case of 
neutron star masses and measurement errors considered 
above. 



IV. DISCUSSION 

In this paper, we have shown that measuring the 
masses and radii of three neutron stars with a 10% uncer- 
tainty can constrain to similar accuracy the pressure of 
neutron-star matter at densities several times the density 
of nuclear saturation. Moreover, astrophysical measure- 
ments with 5% uncertainties can distinguish, to a better 
than 3cr confidence level, not only between equations of 
state that have different underlying physics (e.g., FPS 
and GSl in Figure [5]) but also between ones with very 
similar predicted mass-radius relations (e.g., FPS and 
AP3; see Figure [J). 

Several types of observations in the electromagnetic 
spectrum have already led to a number of mass and ra- 
dius measurements for neutron stars in a wide variety of 
astrophysical systems. Dynamical measurements of the 
masses of radio pulsars in binary systems have uncertain- 
ties as low as 0.1% [1], but do not provide any information 
on the radii of the stars. As such, these observations can 
be used only to set lower bounds on the pressure at the 
three fiducial densities (and especially at p2 and pa) Q. 



Similar constraints, albeit with much larger uncertain- 
ties, can also be set by the dynamical measurements of 
neutron-star masses in X-ray binaries, such as those re- 
ported for Vela X-1 ^\ or Cyg X-2 fg*]. At the same time, 
a secure identification as a black hole of a compact object 
that has a dynamically measured mass of ~ 2Mqi would 
place a stringent upper limit on the same parameters. 

Measurements of both the masses and radii of neutron 
stars have been reported for weakly magnetic neutron 
stars in low-mass X-ray binary systems. In a number of 
cases, measurements were performed for globular cluster 
X-ray transients in quiescence, when the accretion prac- 
tically ceases and the surface emission of the cooling neu- 
tron star can be directly observed [l3, [HI • These obser- 
vations constrain the apparent radii of neutron stars and, 
therefore, lead to measurements of their masses and true 
radii with significant degeneracies. The degeneracies can 
be overcome, however, if the masses of the neutron stars 
are measured independently with dynamical observations 
of the binary systems in which they reside. Current ob- 
servations with XMM-Newton led to ~ 15% uncertainties 
in the measurements of the apparent radii jll| . Using a 
telescope with an effective area of 2 m^, such as the one 
proposed for the International X-ray Observatory, would 
reduce this uncertainty to ^ 5%. 

A different set of measurements of both the masses 
and radii of weakly magnetic neutron stars were recently 
reported that utilized observations of their thermal emis- 
sion during thermonuclear bursts 0, [l2| • The degenera- 
cies between the inferred masses and radii were overcome 
in this case because up to three distinct observables that 
have different dependencies on the masses and radii were 
used. In this case, the measurements were performed 
with the Rossi X-ray Timing Explorer and resulted in 
uncertainties of ~ 10%. Reducing the uncertainties to 
5% requires an increase of the effective area of the tele- 
scope by a factor of ~ 4, i.e., an effective area of 3.2 m^. 
This is well within the design specifications of the Ad- 
vanced X-ray Timing Array (AXTAR). 

The masses and radii of neutron stars can in principle 
be constrained also by observing and modeling in detail 
the fiux oscillations at the spin frequency of the neutron 
stars during the spreading of thermonuclear flashes [l^l . 
This is not currently possible because it requires count 
rates significant larger than those achieved by the Rossi 
Explorer even for the brightest sources. Such measure- 
ments, however, may be achieved with future observa- 
tions. 

Dynamical measurements of the masses and moments 
of inertia of at least three neutron stars may also lead to 
very similar constraints as those discussed in this work 0, 
Indeed, the moments of intcrtia of neutron stars of 
known masses encode not only information about their 
radii but of their density profiles, as well. As such, they 
are sensitive to the parameters of the equation of state 
at low densities and will lead to improved constraints on 
the parameters Pi and po- For the case of the double 
pulsar, a measurement of the moment of inertia with an 
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uncertainty of order 10% will probably be achieved by 
2020 15]. 

Finally, observations of the end stages of the coales- 
cence of double neutron stars with gravitational wave 
observatories such as LIGO, VIRGO, and GEO-600 will 
open new avenues for measuring the masses and radii of 
neutron stars and, hence, for determining the properties 
of the equation of state of ultradense matter [16|] . 
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